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Abstract 

Two different sets of collective coordinate equations for solitary so- 
lutions of Nonlinear Klein-Gordon (NKG) model are introduced. The 
collective coordinate equations are derived using different approaches 
for adding the inhomogeneities as exrernal potentials to the soliton 
equation of motion. Interaction of the NKG field with a local inho- 
mogeneity like a delta function potential wall and also delta function 
potential well is investigated using the presented collective coordinate 
equations and the results of two different models are compared. Most 
of the characters of the interaction are derived analytically. Analyt- 
ical results are also compared with the results of numerical simulations. 
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1 Introduction 



Solitons are localized waves that have a nonzero energy density in a finite 
region of space and exponentially goes to zero as one moves away from 
this region. They appear in nonlinear classical field theories as stable and 
particle-like objects with finite mass and explicit structures. Therefore find- 
ing suitable methods for studying the soliton as a point particle help us 
to find a better perspective of the soliton behaviour. On the other hand, 
comparing the results of such these models with the results of direct nu- 
merical simulations determines the differences between solitons as point-like 
particles and real solitons. This topic is an interesting subject in nonlinear 
field theories pQ. Solitons appear in a nonlinear medium with a fine tuning 
between Nonlinear and dispersive effects. This means that they may be dis- 
appeared in the lack of this precise balance in the medium. It is clear that 
a real medium contains disorders and impurities. Therefore stability and 
propagation of solitons in such these media are of great interest because of 
its applications and also theoretical interests. In order to understanding the 
behaviour of nonlinear excitations in a disordered system, it is important to 
investigate the interaction of solitons with impurities. 

Recently, some non-classical behaviours have been reported for solitons dur- 
ing the scattering from external potentials [2]. These potentials generally 
come from medium defects or impurities. Scattering of solitons of integrable 
systems from the potentials have been studied before [3J; but such this in- 
vestigation for non integrable systems has not been reported yet. Therefore 
it is interesting to examine the methods of adding the potential to the NKG 
model as a non-integrable model and compare the outcomes with the re- 
sults of integrable systems. These are strong motivations for investigation 
the scattering of NKG solitons from defects. 

External potentials can be added to the equation of motion using differ- 
ent methods. One way is that adding an external potential to the equation 
of motion as perturbative terms [2j [3J . These effects also can be taken into 
account by making some parameters of the equation of motion to be func- 
tion of space or time [U [5] . Another way is adding an external potential to 
the field through the metric of background space-time [3 E] • This method 
can be used for models that their lagrangians are Lorentz invariant, such 
as Sine-Gordon model, </> 4 theory, CP N model, NKG models and etc. In 
this paper we will focus on the behaviour of solitons of NKG and try to 
investigate the interaction of NKG solitons with defects by introducing two 
different analytical models. 

Different types of the NKG model are important non-integrable models 
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which appear in some branches of science. These equations can be used 
for describing the particle dynamics in quantum field theory. Some of the 
other examples of NKG applications are: Discrete gap breathers in a di- 
atomic chain [9] , Dichotomous collective proton dynamics in ice [10] , propa- 
gation and stability of relaxation modes in the Landau-Ginzburg model with 
dissipation [11] and Pion form factor |12j . Recently Wazwaz has proposed 
several localized solutions for NKG equations using "Tanh" method |13j . 
Solitons present different trajectories during the interaction with potentials. 
They can pass through or becomes trapped inside the potential after the 
interaction. This behaviour is very sensitive to the values of potential pa- 
rameters in the model as well as to the initial conditions of scattered soliton. 
Most of the researches are in base of numerical studies in nonlinear field theo- 
ries. Because such these systems are generally non-integrable. The collective 
coordinate approach helps us to find analytical equations for the evolution 
of localized solutions, clearly if one can construct such suitable variables. 
We will present two sets of collective coordinate variables extracted from 
different hypotheses [H]. These help us to talk about the validity of their 
results and predictions. 

Therefore two models for the NKG field in an space dependent potential is 
presented in section 2. Both analytical models are introduced and will be 
solved in section 3. The results of two analytical models are compared in 
section 4 for potential barrier and potential well systems. In section 5 we 
will compare our analytical results with direct numerical solutions of the 
equations. Some conclusion and remarks will be presented in the section 6. 

2 Two analytical models for NKG soliton-potential 

system 

Model 1. Lagrangian of the NKG model in (1+1) dimensions is defined as 



C = -d^<j> - U (0) 



(2.1) 



where U {4>) is potential of the field which is defined by 




(2.2) 



The equation of the motion from the Lagrangian (2.1) is 




(2.3) 
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Equation (2.3) has a soliton solution which can be written as [13 

1/2 



4>(x,t) 



, ( prX-Xit) 

1 + tank VA — = 



(2.4) 



where X{t) = xq — Xt . xq and X are soliton initial position and its ve- 
locity respectively. It is a kink- like solution as figure 1 shows. We want to 




Figure 1: Kink-like solution of the NKG described by equation (2.3) for 
xq = and u = 0.5 at t = 0. 



investigate the behaviour of kink solution (2.4) during the interaction with 



an external potential. The external potential V(x) can be added to the La- 
grangian by taking a suitable space dependant function for the parameter A 
as A(x) = 1 + V(x) [HI [151 [16]. I n this approach we have added the potential 
directly to the equation through an extra term. 



By inserting the solution (2.4) in the Lagrangian (2.1) and using adiabatic 
approximation [2] we have 



4 



X 2 - lj sech 4 (x - X) 
8(l + tanh {x - X)) 



\{x) 



1 + tank (x - X) - (1 + tank (x - X)) 2 + 



[l + tank (x- X)f 



(2.5) 



Model 2. The general form of the action in an arbitrary metric is 



S 



C{4),d^)^gd n xdt 



(2.6) 



where "g" is the determinant of the metric g^ v (x). Energy density of the 
system can be found by varying both the field and the metric [6|. For the 
Lagrangian of the form (2.1 ) the equation of motion becomes [U Q3] 



1 



-gd„<t>d»<t> + d^V^g) + 



dU(4>) 







(2.7) 



-g - ■ ■ d(j) 

Space dependent potential can be added to the Lagrangian of the system by 
introducing a suitable nontrivial metric for the back ground space-time 
I14j . In other words, the metric carries the information about the potential. 
In the presence of a weak potential V(x) the suitable metric is [HI El E! 

1 + V(x) 
-1 



(2.8) 



By inserting the solution (2.4) in the Lagrangian (2.1) with the potential 
(2.2) of NKG and using the metric (2.8), with adiabatic approximation [21 E] 
we have 



C = y/l + V(x) ((1 + V{x)) X 2 -l)^ 



sech 4 ((x - X)) 
(1 + tank (x - X)) 



V 1 + W) 



1 + tanh ( x _ x) - (1 + tank (x - X)f + ( X + tanhjx - X)f 

(2.9) 



For the weak potential V(x) (2.9) becomes 

V(x) 



C 



1 + 



3V(x) 



X" 



1 + 



sech 4 {(x - X)) 
■ (1 + tank (x - X)) 



1 + 



V(x) 



1 + tank (x - X) - (1 + tank (x - X)) 2 + i 1 + tanh ( x ~ X )f 

(2.10) 
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3 Collective coordinate for two models 



The Lagrangian density of the soliton is described by (2.5) in the model 



1 and (2.10) in the model 2. These two equations are different in kinetic 



and also potential terms. We will compare them later. The soliton internal 
structure can be ommited by integrating the Lagrangian density (or Hamil- 
tonian density) respect to the variable x. The integrated Lagrangian calls 
collective Lagrangian. After the integration, the soliton spears as a point- 
like particle, however the effects of its extended nature is reflected to the 
kinetic and also potential parts of the collective Lagrangian. The dynamics 
of the point-like particle can be described by equations which are derived 
using collective Lagrangian. It is interesting to compare the results of col- 
lective equations with direct numerical simulation of the main Lagrangian 
density. Let us derive the collective Lagrangian and the point-like particle 
equation of motion in two models. 

Model 1. By integrating Lagrangian (2.5) over the variable x , X(t) re- 
mains as a collective coordinate. If we take the potential V(x) = e8(x) 



Collective lagrangian is derived from (2.5) as 



L 



-X 2 + - 



(tanh (X) - 1) + (1 - tank (X)Y 



[l -tank (X)Y 



(3.11) 



where Mq = g is the rest mass of the soliton in this model. The effective 



potential comes from equation (3.11) which is given by 



tank (X) - 1 + (1 - tank (X)f 



{I -tank (X)Y 



+ 



1 



2 

(3-12) 

NKG solitons have a special and interesting situation which has not been 
observed in other fields. If we plot the effective potential as a function of 
collective position (X), we will find that it has a spatial shift respect to the 
origin. Figure 2(a) shows the effective potential as a function of position (X) 
for e = 2. Simulations also confirm this point for the NKG model. Figure 
2(b) presents the shape of the potential barrier as seen by the soliton which 
is plotted using numerical calculation. It is an special characteristic for the 
NKG field theory. The source of this repugnance with the other fields needs 
more reconsideration. 



The equation of motion for the variable X(t) is derived from (3.11) as 
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X 



-2 

(b) 



o 

position 



2 



Figure 2: (a) Effective potentiai as a function of position with e = 2. (b) 
Potential barrier as seen by the soliton in the model 1 using numerical 
simulations. 



-X - -sech 2 (X) 
2 2 v ; 



1-2(1 -tank (X)) + 



3(1 -tanh (X)f 



(3.13) 



We can define a collective force on the soliton if we look at the above equation 
as F = MX, where M is the rest mass of the soliton. Therefore we have 



F = -sech 2 (X) 



3 11 

-tank 2 (X) + -tank (X) - - 

4 v ' 2 4 



(3.14) 



The above equation shows that the peak of the soliton moves under the 
influence of a complicated force which is a function of external potential 
and soliton position. Suppose that the soliton moves toward a potential 
barrier. Its velocity reduces because of affecting a repulsive force. After 
that soliton velocity increases when it goes away. Figures 3(a) and 3(b) 
show the force exerted by the potential well and barrier on the soliton for 
e = —4 and e = 4 respectively. This is in agreement with the observed 
behaviour for NKG model as shown in figure 2. These figures also show 
that the center of the f orce i s not located in the origin. 
Fortunately equation (3.13) has an exact solution for X as follows 



X 2 -X 2 = e - [tank 3 (X) + tank 2 (X) - tank (X) - tanh 3 (X ) 



tank 2 (X ) + tanh (X )] 
(3.15) 
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Figure 3: (a) The force on the soliton by a potential well with e = —4. (b) 
The force on the soliton by a barrier with e = 4. 



where Xq and Xq are soliton initial position and its initial velocity respec- 
tively. Some of the physical features of soliton-potential system can be 



found using equation (3.15 ). Collective energy is obtainable from Lagrangian 



(3.11 ) as follow 



E 



:X l 



{tank (X) - 1) + (1 - tank (X)Y 



{I -tank (X)Y 



+ 



1 



2 

■ (3 - 16) 

It is the energy of a particle with the mass of Mq = ^ and velocity X which 
is moved unde r the influence of exter nal effective potential. By substituting 
X from (3.15) in the equation (3.16), one can show that the energy of the 



system is a function of soliton initial conditions Xq and Xq and therefore it 
is conserved. 

Model 2. The above results can be calculated using the model 2. For 



the Lagrangian (2.10) X(t) remains as a collective coordinate if we inte- 
grate (2.10) over the variable x. If we take the potential V(x) = e5(x), the 



collective Lagrangian becomes 



L 



1 3esech 4 (X) 
4 + 16(1 



X 1 



+ - 



(tank (X) - 1) + (1 - tank {X)) z 



tank (X)) 
(l-tanh(X)f 



sea 



h*(X) 



4(1 



tank (X)) 
(3.17) " 



S 



The equation of motion for the variable X(t) is derived from (3.17) 



+ -sech 2 (X) 
4 



1 ^ 3esec/i 4 (X) 

2 + 8 {1- tank (X)) 
1 



X 



-tanK 1 (X) - -tank (X) + 1 



3X 2 + 2) = (3.18) 



The above equation describes the soliton trajectory which moves under the 
influence of a collective force. The collective force is a function of soliton 
position and also its velocity. X{t) can be calculated as a function of X(t) 



by integrating (3.18) as follows 



3X 2 + 2 
3X 2 + 2 



+ 



+ 



3esech 4 (Xo) 
8(l-tanh(X )) 

3esech 4 (X)~ 
8(l-tanh(X)) 



(3.19) 



Where Xq and Xq are soliton initial position and its initial velocity respec- 
tively. The energy of the soliton in the presence of the potential V(x) = e5{x) 
using model 2 becomes 



E 



1 3esec/i 4 (X) 
4 + 16(1 



X< 



{tank (X) - 1) + (1 - tank {X)Y 



tank (X))_ 
{I -tank {X)f 



sech 4 (X) 



4(1 



tank (X)) 
(3.20) 



+ 



Equation (3.20) shows that the rest mass is a function of soliton position in 



the model 2. This is the reason of some differences between two models which 
will be discussed in the next section. Some features of the soliton-potential 



dynamics are also studied using equations (3.19) and (3.20) analytically in 
the section 4 too. 



4 Comparing the models 

Potential barrier. There are two different trajectories for a soliton during 
the interaction with an effective potential barrier which are depend on the 
soliton initial conditions. A soliton with a low velocity reflects back from the 
barrier and a high velocity soliton climbs over the barrier and passes over 
it. So these two situations can be separated by a critical velocity. The total 
energy of the soliton-potential is conserved in both two models, as mentioned 
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before. Therefore we can find the critical velocity with a simple analysis from 
the energy of th e solit on w ithout any needed numerical simulations. Both 

\Xn + if when the 



equations ( 3. 161) and (3.20) reduce to E(X = oo 



soliton is far from the center of the delta-like potential which is located at 

} and velocity 



the origin. It is the energy of a particle with the mass Mq 
of Xq . The energy of a soliton in the origin (X = 0) comes from (3.16) and 
( 3.20 ) for two models: Ei(X = 0) = \Xq + 
E 2 (X = 0)- n 



\ for the model 1 and 



>4 ^ 16 



) Xq + 5 for the model 2. The minimum energy of 
soliton in this position is E = | + ^ for two models. On the other hand, a 
soliton which comes from the infinity with initial velocity v c has the energy 



of E 



+ 



2 4 



So it is easy to calculate the critical velocity of soliton by 



comparing the energy of the soliton at the origin and infinity. The critical 
velocity is obtained as v c = a/J using both two models. The same result is 
derived by substituting X = , Xq = v c , Xq = oo and X = in (3.15) and 



(3.19). 



Note that the critical velocity of the soliton depends on its initial position 
as well as its initial velocity. For a soliton which is located at some position 
like Xq (which is not necessary infinity) the critical velocity will not be 
So a soliton in the initial position Xq with initial velocity of Xq 



2" 



has the critical initial velocity if its velocity becomes zero at the top of the 
barrier X = 0. Consider a soliton wit h init ial co nditio ns of Xo and Xq. If 
we set X = and X = in equations (3.15) and (3.19) then v c = Xq. Thus 
for the model 1 we have 



(tanh (Xq) - tanh 2 (X ) - tanh 3 (X )) (4.21) 



But the critical velocity in the model 2 becomes 



'2e (1 - tanh(Xo) - sech 4 (X )) 
4 - 4tanh(X ) + 3esech 4 (X ) 



(4.22) 



Figure 4 shows critical velocity as a function of potential strength for Xq = 
— 1 in two models. This figure shows that the model 2 predicts smaller 
critical velocity respect to the model 1 because of differences between rest 
masses of these models. It is clear that a soliton with great rest mass needs 
smaller velocity to reach the potential peak. Difference between the rest 



(3.20) as 



mass of the soliton in the models 1 and 2 can be calculated using (3.16) and 

3esech A (X) 



AE k 



inetic 



16 (1 — tanh(X)) 



X' 



(4.23) 
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The above equation shows that two models predict equal rest mass at the 
infinity. But the difference between calculated rest mass in two models 
increases when the strength of the potential increases. Figure 4 shows this 
phenomenon explicitly. It is interesting to depict the critical velocity as a 




0.2- 



1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

Barrier height 



Figure 4: Critical velocity as a function of barrier height in both models for 
initial position Xq = — 1. 

function of initial position. The critical velocity to pass over the potential 
has been demonstrated as a function of the initial position in the figure 5 
for two models with e = 0.5. This figure shows there is a good agreement 
between these two models. For a soliton comes from infinity models are 
equal to each other as one can find from the figure 5. But there are some 
differences for a soliton moves from an initial posotion closer to the center 
of the potential. Figure 5 also demontrates that soliton needs lower initial 
velocity to pass over the barrier if it is closer to the center of the potential. 

Soliton-well system. Let us consider a soliton which moves toward a 
frictionless potential well. This situation is very interesting because of some 
differences between point particle and a soliton in the potential well. A 
point particle falls in the well with an increasing velocity and reaches the 
bottom of the well with its maximum speed. After that, it will climb the 
well with a decreasing velocity and finally passes through the well. Its final 
velocity after the interaction is equal to its initial speed. 



Potential well can be obtained by replacing e with — e in the equation (3.15). 
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Figure 5: Critical velocity as a function of initial position in both models 
for e = 0.5. 

The solution for the system in the model 1 is 

X 2 -X 2 = -| [tanh 3 (X) + tank 2 (X) - tanh (X) - tank 3 (X ) - tank 2 (X ) + tanh (X ) 

(4.24) 

Similary for the model 2 we have 

■ „ 1 3esech 4 (Xo) 

+2 _ 2 _ 8(l-tanh(X )) , , 

QV 2 _LO 1 3esech4(X) 

6J^o -+- Z 2 8 (l-tanh(X)) 

We can define an escape velocity instead of critical velocity for a soliton-well 
system. The escape velocity is the minimum velocity for a soliton which can 
pass through the well. A soliton in an initial position Xq reaches the infinity 
with a zero final velocity if its initial velocity is 



X escap ei = \j- {sech 2 (Xq) - tanh 3 (X Q ) + tanh (X )) (4.26) 



X - / 2esechHX ) 

escape2 ~ y 4 - 4tanh(X ) - 3esech*(X ) { ' 

calculated using the models 1 and 2 respectively. In other words, a soliton 
which is located in the initial position Xq can escape to infinity if its initial 
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velocity Xq is greater than the escape velocity X escape . Figures 6 and 7 show 
escape velocity from a potential well as a function of the well depth (figure 
6) and soliton initial position (figure 7) using two models. A soliton in the 
model 1 needs lower escape velocity in comparson with the model 2 due to 
its bigger rest mass as one can find from the equation (4.23). 
Consider a soliton which moves toward the potential well with an initial 




Figure 6: Escape velocity as a function of the well depth in both models for 
initial position Xq = — 1. 

velocity Xq smaller than the escape velocity X escape . The soliton reaches a 
maximum distance X max from the center of the potential with a zero velocity 
and then come back toward the center of the potential well. Therefore the 
soliton oscillates in the well with the amplitude of X max . The required initial 



velocity to reach X max is found from (4.25) for the model 2 as 
Xq~- 



'2 (1 - tank {Xq)) (4 - 4tanh(AW E ) - 3esech A (X max )) 2 



^ 3(l-tanh(X ma;c ))(4-4tanh(Xo)-3esec/i 4 (Xo)) 3 

(4.28) 

It is clear that the soliton oscillates around the well if its initial velocity is 
lower than the escape velocity. The period of the oscillation can be calcu- 



lated numerically using equation (4.25). 
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Figure 7: Escape velocity as a function of initial position in both models for 
e = 0.5. 

5 Analytical results versus numerical simulation 

It is shown that the general behaviour of a soliton-potential system is almost 
the same in the models 1 and 2. However, we can find some small differences 
between the dynamics of a soliton in different models. It is important to 
compare the results of analytical models with direct numerical solutions. 
Here we will compare the results of analytic model 1 with its numerical 
solution. It is clear that the same comparison can be done for the model 2. 
The soliton equation of motion in the model 1 for a small potential is [15] 



dU 

<fa ~ 4>xx + (1 + <tS{x)) — = (5.29) 

ocp 

Both two models (for a delta-like potential) have a parameter in their equa- 
tion of motion which controls the strength of the external potential. It is 
possible to compare the strength parameters in a specific situation by sim- 
ulation and adjusting parameters to have same results by different models 
for that specific situation. It is expected to find approximately same rela- 
tion between the parameters in other situations. We found that the critical 
velocity for a soliton-barrier system is v c = It is possible to adjust 



the strength parameter e in analytical model with same parameter in (5.29) 
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trough the v c . Numerical simulations using equation (5.29) show the same 



behaviour for critical velocity. An effective potential can be found by inter- 
polation of simulation results on the v c = \ ^P- = \ a + e P . £e , , i s potential 



strength (as an effective parameter) in numerical simulation while e is sim- 
ilar parameter in the model 1. e e ff can be found by fitting the numerical 
results on a theoretical diagram as 



(0.0434 ± 0.01061) + (0.76462 ± 0.02479) e 



(5.30) 



^effective 

Figure 8 shows the result of simulations of equation ( |5.29 ) for the NKG 
model. Our simulations show very good agreement between numerical re- 
sults and theoretical predictions for other features of soliton-potential inter- 
action. 



re 0.3- 



=0.0434+. 76462*c 



analytical 



0.4 0.6 

Barrier height 



— i 

1.0 



Figure 8: Critical velocity as a function of e with results of simulation using 



equation (5.29) and analytic model. 



6 Conclusion and remarks 

Two analytical models for the interaction of the NKG solitons with delta 
function potential have been presented. Models predict a critical velocity 
for the soliton-barrier interaction which is a function of initial conditions 
and the potential identities. For a soliton-well system an escape velocity 
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was introduced instead of the critical velocity. These models are able to 
explain most of the features of the system analytically. We have observed 
that the center of the potential as seen by NKG solitons is a quite different 
from the real position of the potential center. Numerical simulations are 
in agreement with theoretical predictions of the models. Our models fail 
to predict the narrow windows of soliton reflection from the potential well. 
So, it is expected to find a better model with suitable collective coordinate 
method to explain this behaviour. These models can be used for prediction 
of soliton behaviour in the other field theories beside the NKG model. 
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